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a^ ■ ABSTRACT 

jlLi' We derive constraints on the strength of a new intermediate range interaction that couples to 

1^ • baryon number from primordial nucleosynthesis yields. The nucleosysnthesis limits here used 

- ^ , arise from matching observations and predictions of standard and inhomogeneous primordial 

5h I scenarios. We show that the standard nucleosynthesis scenario is more restrictive (as < 0.2) 



when the range of the interaction is greater than about 1 m. We further discuss the implications 
of considering the scalar particle responsible for the new interaction as the main component 
of the dark matter in the galactic halo such that its decay can account for the ionization of 
hydrogen in the interstellar medium and the temperature of Lyman-a clouds. 



1. The existence of new fundamental forces beyond the aheady known four interactions 
is an exciting possibihty that may have profound imphcations in our understanding of the 
physics beyond the standard model. The claim, more than a decade ago, of evidence for an 
intermediate range interaction with sub-gravitational strength [^ has led to a great demand 
of theoretical explanations (see [0 for a review and a complete set of references) and, most 
importantly, has given origin to fresh experiments based on new ideas and to the repetition of 
well known experiments using new state of the art technology. 

In its simplest versions, the putative new interaction or a fifth force would arise from the 
exchange of a light boson coupled to matter with a strength comparable to gravity. There are 
several schemes through which physics at the Planck scale could give origin to such an inter- 
action and yield a Yukawa type modification in the interaction energy between point masses. 
This new contribution to the interaction energy can arise, for instance, from extended super- 
gravity theories after dimensional reduction 0, Q , from the compactification of 5-dimensional 
generalized Kaluza-Klein theories that include gauge interactions at higher dimensions [^ and 
also from string theory. On quite general terms, the interaction energy, V{r), between two 
point masses mi and 1712, can be expressed in terms of the gravitational interaction as 

ir/ \ Goo nil 1712 /. , -r/A5^ n\ 

V{r) = (l + age ' ") , (1) 

r 

where r = |r*2 — n| is the distance between the masses. Goo is the gravitational coupling for 

r —* 00, a^ and A5 are the strength and the range of the new interaction. Of course. Goo has to 

be identified with the Newtonian gravitational constant and the gravitational coupling would 

be dependent on r. Indeed, the force associated with eq. (|ID is given by: 



Fir) = -VV(r) = -^M^ f . (2) 

where 



Gir) = Goo[l + as (1 + r/As) e-'"/"'^] . (3) 



The great interest sparked by the suggestion of existence of a new interaction was the recog- 
nition that the couphng a^ was not an universal constant, but instead a parameter depending 
on the chemical composition of the test masses |5|. This dependence comes about if one as- 
sumes that the new bosonic field couples to the baryon number B = Z + N which is the sum 
of protons and neutrons. Hence the new interaction between masses with baryon numbers Bi 
and B2 can be expressed through a new fundamental constant, /, in the following way: 



V{r) = -f ^^ e-^/^^ , (4) 



such that the constant 0:5 can be written as 

with a = P/Goo^H ^^^ /^i,2 = '"^i,2/"^_H" i^TT'H IS the hydrogen mass). 

Of course, from the above equations it follows that in a Galileo-type experiment an accel- 
eration difference between masses rrii and 1712 would be given by: 



(B' 



'BA _ (B^ 
where F is the field strength of the Earth (which is denoted by 



(6) 



Several experiments (see, for instance, [0] for a list of the most important ones) have been 
performed in order to establish the parameters of a new interaction based on the idea of a 
composition-dependence differential acceleration as described in eq. (|^) and other composition- 
independent effects. The current experimental situation is essentially compatible with predic- 
tions of Newtonian gravity in either composition-independent or composition-dependent exper- 
iments. The bounds on parameters a^ and A5 can be summarised as follows: 

- satellite tests probing ranges about 10^ m < A5 < 10^ m indicate that a^ < 10^^; 

- gravimetric experiments that are sensitive in the range of 10 m < A5 < 10^ m suggest that 
as < 10-^ 

- laboratory experiments deviced to measure deviations from the inverse-square law are sensitive 
essentially to the range 10~^ m < A5 < 1 ?7i and constrain a^ to be smaller than 10~^. 



Interestingly, for A5 < 10~^ m and A5 > 10^^ m, a^ is essentially unconstrained. The former 
range is particularly attractive as new forces with sub-millimetric range seems to be favoured 
from scalar interactions in supersymmetric theories [^ and in the recently proposed theories on 
TeV scale quantum gravity 0. This range also arises from assuming that scalar ||^ or tensor 
interactions associated to the breaking of the Lorentz invariance in string theories account 
for the vacuum energy up to the level fly < 0.5. 

In order to close our summary of the experimental situation we should point out that, as 
discussed in [^, existing experimental data cannot account for certain anomalies such as the 
one claimed to exist in the original Eotvos experiment ||10[ and hence these anomalies remain 



still an open issue. Another quite recent claim concerning the existence of a new interaction, is 
given by the radio metric data from Pioneer 10/11, Galileo and Ulysses spacecraft, indicating 
an anomalous constant acceleration of about 8.5 x 10~® cm s~^ acting on the spacecraft directed 
towards the sun. A new interaction with a^ = —1 x 10^^ and A5 = 200 A.U. = 1.49 x 10^^ m 
seems to account for the anomaly ||lT| (see however ||1^). 



In this work we shall establish limits to the coupling 0:5 and range A5 of the putative new 
interaction using results from nucleosynthesis in the context of the Big-Bang. As we shall see, 
these limits appear to be much less stringent than the ones obtained from laboratory exper- 
iments. Nevertheless, independently from the theoretical setting from which this interaction 
may arise, if the coupling of the putative new interaction is a running coupling constant, the 
limits derived here may turn out to be relevant. Indeed, in this case, the coupling constant 
at the time of primordial nucleosynthesis af""^ can be greater than the bounds arising from 
laboratory experiments, as long as the interaction is a local gauge interaction and its /^-function 
is positive. This is, for instance, the case of the U{1) coupling constant in QED. Assuming 
scale-dependence for the new interaction coupling constant is quite a natural assumption, as 
all gauge coupling constants in the standard model are running. Even the gravitational cou- 
pling constants are, at one-loop level, running couplings in higher derivative theories of gravity 
[|l^] (including of course the Newton's constant). This fact has implications for quite a few 
diversified problems: the problem of the rotation curve of galaxies 0, |T^; the cosmological 
dark matter problem |]T^ , the large scale structure of the Universe |T^ and the cosmic virial 



theorem [1^. Notice that we are concerned only with zero temperature running effects as those 



for couplings of the type gipipip (cf. eq. (|14|)). The finite temperature corrections that are 
proportional to g^ (see eg. |]T9[) become negligible (cf. considerations after eq. (|15|)) in the 
limit T » 1715 (since we shall deal with temperatures T ~ MeV and argue that m^ ~ eV). 

Of course, the limits that we are going to establish can be regarded as independent of any 
considerations concerning the running of couplings and are on their own of relevance as they are 
consistent with laboratory experiments and are obtained from an independent line of reasoning 
than the usual approaches. 



Standard primordial nucleosynthesis scenario (see [^ for a review) allows obtaining bounds 
on the variation of the effective gravitational coupling when confronted with the measurements 
of abundance of light elements in the Universe. Such measurements are in agreement with the 
standard nucleosynthesis scenario leaving still some room for variations in the effective number 
of neutrinos, the baryon fraction of the universe and also in the value of the gravitational con- 
stant. It is well known that the predicted mass fraction of primordial "^He can be parametrised. 



in theories with a variable gravitational coupling, in the following way |20, 21 



rp = 0.228 + 0.010 In r/io + 0.327 log ^ , (7) 

where rjio denotes the baryon to photon ratio in units of 10~^° and ^ is the ratio of the Hubble 
parameter at the time of nucleosynthesis and its present value which is itself proportional to 
the square root of the gravitational constant. In the fit eq. (^ it is assumed that the effective 
number of light neutrinos is N^ = 3 and that the neutron lifetime is r„ = 887.0 ± 2.0 seconds 



22 . 



A range for the values of the effective gravitational coupling that are compatible with the 
observations of the primordial D, '^He, ^He and ^Li abundances can be obtained running the 



nucleosynthesis codes for different values of G ||2^ . It turns out that the permissible range is 



rather large {AG/G = 0.2 at the 1 a level), given the large statistical and systematic errors of 
the observations. We shall use this result to constrain the parameters of a new intermediate 
range interaction assuming it to be sensitive to baryon number. The maximum range that can 



be tested corresponds to the physical horizon distance at the time of primordial nucleosynthesis. 
At that time, the horizon distance grows from a few light-seconds to a few light-minutes, and 
therefore is much smaller than a few microparsecs: r < 3x 10^° m. Using eq. (13) and considering 



the 1 a level result of Ref . p3[ , we obtain: 



AG{r 



^ = as (1 + r/As) e""/^^ < 0.2 



a. 

which we shall use in the next sections to impose constraints on the parameters of the new 
interaction. 

Of course, as a light-second is about the distance to the Moon, constraints on a variation of 



G can, at this scale, be set from lunar laser ranging, AG/G < 0.6 p4|. It is worth mentioning 



that similar reasoning has been used to impose constraints on the variation of the gravitational 



coupling due to a possible dependence on scale P 



2. Let us turn to the discussion of the bounds that can be imposed from the standard 
primordial nucleosynthesis scenario. First of all we need to estimate a typical distance r between 
particles interacting through this new interaction. Only then can we infer on the range and 
strength of the fifth force. 

The period of interest is that immediately after the weak-interaction decoupling, for which 
temperatures are of the order E' ~ 1 MeV and the total baryon density is pb ~ 1.21 x 
10^^ g cm^^ ||2^. We shall assume that, in this epoch, the baryons are in the form of protons 



and neutrons only: ps = Pp + Pn- 

One can argue that the typical distance between the interacting particles should be smaller 
or of the same order as their mean free paths A„ and Xp. In order to estimate A„ and Xp 
we need to calculate the relevant stopping cross sections and the densities of the stopping 
particles. Within a standard scenario the number of neutrons is slightly smaller than the 
number of protons due to the mass difference and neutron decay: 

n^ Am c^ t 1 

— = e fcT e ^" ~ - , (9) 

Up e 

although, at this temperature t << Tn, and hence the neutron decay is not so important. 



Using eq. (|) and the charge neutrahty condition, we arrive at the following estimates: Pp = 
0.9 X 10~^ g cm~^, pn = 0.3 x 10"^ g cm~^, and pe = 0.5 x 10"^ g cm~^. 

In order to estimate the relevant cross sections we should keep in mind that the proton is 
essentially scattered by the electrons, whereas for the neutron, both proton and electron scat- 
tering processes should in principle be taken into account. For these cases we have reproduced 
the estimates given in [EB]. 



Let us first consider the proton-electron scattering. Here we use the expression from p6| : 

At -E ~ 1 MeV, the Debye shielding of the proton by the electronic cloud inhibits scattering 
below 9q ^ 0.77". Evaluating the integral in eq. ([T0|) we obtain ape — 1.5 x 10~^^ cm^. 

The nucleon-nucleon cross section at these energies is very well known given the large amount 
of data available for the nucleon-nucleon phase shifts. Typically cr„p is written as the sum of a 
singlet and a triplet contribution: 

"^^ ' {asky + {l-0.5rsask^y{atky + {l-0.5rtatk^y' ^ ' 

where the parameters for scattering lengths and scattering radii are obtained from fits to the 
data. We have used a^ = —23.71 fm; r^ = 2.73 fm for the singlet component and at = 5.432 



fm; rt = 1.749 fm for the triplet component |2^. This yields the following value for the 
nucleon-nucleon cross section: anp c^ 1.9 x 10~^^ cm'^. 

The neutrons interact with the electrons through their magnetic moment. At these energies 
(E' ~ 1 MeV) the cross section is given by ||26|| : 



a^^ = 37T (^^^I^f = 8 X 10-31 cm^ , (12) 

where K^ag = —1.91 is the anomalous magnetic moment of the neutron in nuclear magnetons. 

The mean free path is defined as: A = ^^ where m and p are the mass and density of 
the stopping particles and a is the relevant stopping cross section. Taking into account our 
estimates for the cross sections, we easily conclude that, in this epoch, both neutrons and 
protons have mean free paths of the same order of magnitude: Ap ~ A„ ~ A = 1 m. 



In Figure (|T]) we present a contour plot of -pr^ as a function of the interaction parameters 
(asjAs) by taking the typical distance between nucleons to be r ~ A = 1 ?ti. The black, 
dark grey and light grey regions correspond to ^^' within the ranges [0,0.1], [0.1,0.2], and 
[0.2, 0.3] respectively. One can immediately conclude from the contour shapes that, whichever 
the range for the interaction, the maximum restriction that can be extracted for the coupling 
constant as from nucleosysnthesis corresponds exactly to the limit given for q ■ It is worth 
underlining that the condition on a^ is valid for the primordial epoch (given that primordial 
nucleosysnthesis was used). If as is assumed to be a running coupling constant, the hmits 
here derived for the strength of this new interaction are quite useful, given that laboratory 
experiments can never impose constraints on af*^**". 

Evidently af^*™" < 0.2 is far above the typical values permitted by experiments (as" °^ < 
10~^). In spite of that, even if the coupling constant is not running, it is interesting to look 
at the extremes where there are no experimental bounds. If this new interaction is very short 
range As < 10""^ m, primordial nucleosysnthesis does not constrain the coupling constant at 
all. However if this is a very long range interaction (As > 10^^ m), we find that af"™ < 0.2 in 
order to satisfy eq. (||). 

So far we have approximated the typical distance between the interacting particles by the 
mean free paths r ~ A„ ~ Xp. Alternatively, a very trivial estimate for the typical distance 
between nucleons could be extracted directly from the total baryon density normalisation con- 
dition: 

« = i^y^' (13) 

PbNa 
where Na is the Avogadro's number. Evaluating eq. (|13|) yields a ~ 10^^ m. We find that now 

the limit af*^*™ < 0.2 holds even if the interaction is short range a < As < 10"'^ m. Note again 

that there are no experimental limits on as° "^ in this region. 

Ultimately we shall consider the reasoning of the previous section to obtain bounds on the 
parameters of a new intermediate range interaction in the context of inhomogeneous nucle- 
osynthesis scenario. In inhomogeneous Big Bang models one assumes that there are non-linear 
perturbations of the baryon density producing large proton deficient regions where the produc- 
tion of ^He is hindered [^ . The main advantage of these models is that they allow for a larger 
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deuterium abundance, which is one of the directions suggested by observational data |28 . 



In this picture most neutrons drift into neutron rich regions, and therefore we consider 
that it is the electron-neutron scattering process (rather than the proton scattering) which 
determines the neutron mean free path. For the sake of our argument, we shall assume that the 
electrons remain homogeneously distributed. Using the values of cr^e, ape and pe given in the 
previous section, we estimate a neutron mean free path of X^^^"'"^ ~ 10^ m whereas the proton 
mean free path remains essentially the same Ap ~ 1 m,. If we now take the typical distance for 
interacting particles to be r = A^"'^°™ the restrictions imposed on the primordial value of the 
coupling constant are less severe. This aspect is manifest in the contour plot shown in Figure 

(I)- 

We can then conclude that for A5 > 1 m, the standard nucleosynthesis scenario is more 
effective in constraining a^ then the inhomogeneous models. 

3. An independent estimate of the parameters of the new interaction can be obtained 
assuming, for instance, that the light boson carrier of the new interaction accounts for the 
main contribution of the dark matter in the galactic halo, as well as being responsible for 
the ionization of interstellar hydrogen. Even though this might seem highly speculative, the 
clumping of this bosonic field around and inside stars derives naturally from the coupling of this 
field to ordinary matter (cf. eq. ([I^) below). As discussed in Ref. p^ this implies that masses 
and couplings of ordinary particles depends on the features of the bosonic field interaction, 
meaning in turn, that the fundamental coupling constants are altered by the nearby density of 
matter. Estimates of the effects of this dependence on the cooling of neutron stars, the neutrino 
burst in the supernova SN1987A and the period of the remnant pulsar, imply bounds for 0:5 and 
A5 that are much less stringent than the ones emerging from Eotvos-type experiments and from 
satellite measurements mentioned in the introduction. However, as we shall see, interesting 
limits do arise if one assumes that the boson responsible for the new interaction decays into 
photons and accounts for the dark matter contained in the halo of our galaxy. 

Denoting the light scalar field responsible for the new interaction by 0, we assume its 
coupling to nucleons and photons is of the following form: 



(05) (05) 

where (^s) is a large scale associated to the new interaction. Of course, it is this scale that 
establishes the likelihood of creating this field in particle physics accelerators. This interaction 
yields, at tree level, a modification to the Newtonian potential such as that in eq. ([1|) where 

and Mp = G^'^ is the Planck mass. We see that the existing bounds on as imply that 
(05) ~ Mp. This agrees with what had already been mentioned in the introduction: the new 
interaction must arise from physics close to the Planck scale. For simplicity we shall assume 
that (05) = Mp. Following from our previous conclusion that af"*" < 0.2 for A5 > 1 ?7i, we 
obtain c„ < 1.58. 

In order to extract further information on the new interaction, we demand that the scalar 
field, 0, decays into photons that are energetic enough to account for the observed ionization 



of interstellar hydrogen, the temperature of Lyman-a clouds |30|, and the anomaly in the 
abundance of He I in the three high-redshift Lyman-limit systems of the quasar HS 1700 + 6416 
[^] . This implies that m^ > 27.2 eV and hence that A5 < 7 x 10^^ m. This is consistent 
with the bounds on 0:5 obtained from supposing that the distance between nucleons during 
nucleosynthesis is given by eq. (]T3|). 

Assuming the scalar particle is stable enough to account for all the dark matter in the 
galactic halo, then the density of scalar particles is given by p^ = m^n^ = ph, where n^ is the 
scalar particle number density and ph = 2 — 13 x 10~^^ g an~^ = 0.1 — 0.7 GeV cm^^ |2^ is 
the galactic halo density. This hypothesis yields that: 

3.67 X 10^ cm'^ < n^ < 2.57 x 10^ cm"^ . (16) 

Notice that if we had chosen to account for the observed cosmological energy density, 
we would have had to replace ph by (0.2 — 0.3)pc, where pc = 1.88 x 10^^^ /iq g cm~^ = 



10 



1.05 X 10~^ /iq GeV cm~^ is the critical density and Hq is the Hubble parameter in units of 
100 km s~^ Mpc~^ (observationally 0.5 < Hq < 0.8). Then the number density of scalar 
particles would be smaller than the estimate eq. (|1^) by about four orders of magnitude. 

To further verify our assumptions we compute the rate of decay of the scalar particle. Since 
it can only decay into photons, its decay width can be obtained from the last term in eq. ( p!4D . 
A straightforward calculation reveals that: 



For ms = 27.2 eV, we have tg = T^^ = 3.0 x 10=^^ s » tu ^ H^^ = h^^ 9.78 Gyr. Thus 
the scalar particle responsible for the new interaction is comfortably stable to be a good dark 
matter candidate. 

4. In this paper we have discussed bounds on the parameters of a new interaction that 
couples with baryon number, arising from primordial nucleosynthesis. We have considered 
the standard nucleosynthesis scenario and shown that 0:5 < 0.2 for A5 > 1 m (see Figure 1). 



The inhomogeneous nucleosynthesis scenario of Ref. ||2^ allows one to conclude the same for 
A5 > 5 X 10^ m (Figure 2). These limits do not require the assumption of a running coupling 
constant and even though consistent, they are much less stringent than the ones obtained from 
laboratory experiments and satellite tests. As discussed in the text, they may be relevant 
for a new interaction whose coupling constant is scale-dependent and is not asymptotically 
free. We have also derived bounds on the dumpiness of the scalar particle assuming it is the 
main contributor to the dark matter in the galactic halo, that is 3.67 x 10^ cm"^ < n^ < 
2.57 X 10"^ cm~^ for A5 < 7 x 10"^ m. This corresponds to m^ > 27.2 eV implying that the 
scalar particle decay can account for the ionization of hydrogen observed in the interstellar 
medium and the temperature of Lyman-a clouds. 
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Figure 1: Limits imposed on the new force 
parameters based on the nucleosysnthesis con- 
straints for -pr^ within the ranges [0,0.1], 
[0.1,0.2], and [0.2,0.3] respectively. 
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Figure 2: Limits imposed on the new force 
parameters based on the nucleosysnthesis con- 
straints and the inhomogeneous primordial nu- 
cleosysnthesis model for „ ^^' within the ranges 
[0,0.1], [0.1,0.2], and [0.2,0.3] respectively. 
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